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Abstract 

On the level of an effective quark theory, we define confinement by the absence 
of quark anti-quark thresholds in correlation functions. We then propose a 
confining Nambu-Jona-Lasinio-type model. The confinement is implemented 
in analogy to Anderson localization in condensed matter systems. We study 
the model's phase structure as well as its behavior under extreme conditions, 
i.e. high temperature and/or high density. 
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1 Introduction 



The description of hadron physics starting from quantum chromodynamics 
(QCD) - the theory of strong interactions - is one of the most challenging problems 
in medium energy physics of today. The difficulty in the description of the low 
energy sector of QCD lies in the large effective coupling constant, which precludes 
a perturbative treatment. 

Various non-perturbative methods shed some light on its ground state proper- 
ties: the operator product expansion [|l| improves the perturbative expansion by 
the inclusion of effects from vacuum condensates in Green's functions, QCD sum 
rules relate these condensates to hadron masses, variational methods attempt 
to evolve a qualitative picture of the gluonic vacuum, the semi-classical evolution of 
the QCD functional integral, leading to instanton physics H|, and the reformulation 
of QCD in terms of the field strength |5], |], |7|, |8| provides some insight into the 
vacuum properties of the quark sector. These analytical approaches can be con- 
trasted to numerical investigations of lattice QCD || |10[. The latter approach is 
not bounded by any approximation, but is restricted by the capacities of computers. 

One of the most striking features of low energy QCD is the absence of quarks 
in the asymptotic states. This quark confinement is explained in lattice QCD by a 
linear rising potential between two static quarks. This behavior is also confirmed 



by considerations within the i/iV c -expansion, N c being the number of colors |TT 
A natural explanation of the linear confining potential is provided by the dual su- 
perconductor picture [|12|]. As pointed out by 't Hooft, in a certain gauge, the non 



abelian Yang-Mills theories possess monopole configurations |13[ . If these monopoles 



condense, a dual Meissner effect occurs, expelling electric field strength out of the 
vacuum. This implies that the electric field between two static quarks is squeezed 
into a flux tube, which subsequently gives rise to the linear confining potential. A re- 
cent SUSY model of Seiberg and Witten suggests that monopole condensation could 
be the mechanism of confinement in certain Yang- Mills theories ]H , |T5| . There have 



been suggestions that the same phenomenon occurs in the Yang-Mills sector of QCD 
proper In this development, it is found that the confinement of the quarks is 



intimately related to the spontaneous breakdown of chiral symmetry [|Tj 

Despite the recent progress in understanding the ground state properties of QCD, 
the description of hadron properties is still feasible only with effective models. These 
models include aspects of QCD by incorporating its symmetries. The most impor- 
tant symmetry constraining the variety of hadron models is chiral symmetry. For 
QCD (N c = 3) and for zero current masses, the quark sector is invariant under 
global SU(Nf)i x SU(Nf)n transformations of the left- and right-handed quarks, 
where Nf is the number of quark flavors. The vaccum breaks this symmetry down 
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to the diagonal SU(N f )$ 

Most effective theories of mesons and nucleons contain many parameters restrict- 
ing the predictive power of these models. It is therefore desirable to "derive" the 
effective hadron theory from an underlying effective quark model in order to obtain 
constraints on the parameter range. One of these quark models is the model of 
Nambu and Jona-Lasinio (NJL) |]I8|| . The NJL model is one of the most econom- 
ical models that possess the essence of dynamical symmetry breaking that is the 
hall-mark of modern field theories and has enjoyed an impressive phenomenological 



success in hadron physics [19|. Now the NJL-model reflects the principal low-energy 
QCD symmetry properties of the quark sector, but does not include the confinement 
of quarks. This implies, in particular, that the mesons of this model - quark-anti- 
quark bound states - can decay into free quark-anti-quark pairs, if such a process is 
allowed by kinematics. This manifestly unphysical threshold puts severe constraints 



on the applicability of the NJL-model [21 



In this paper, we propose an NJL-type model that possesses the confinement 
property in the sense that quark-anti-quark thresholds are absent in (mesonic) 
Green's functions. Our reasoning, albeit exploratory, will be guided by a close anal- 
ogy to a phenomenon in condensed-matter physics known as Anderson localization 
of electrons [p2]1 . In Anderson localization, it is observed that freely moving electrons 
get localized when the strength of a random distributed potential stemming from 
the impurities of the conducting solid exceed a certain strength. We shall argue in 
analogy that the quarks feel randomly distributed background fields generated by 
the gluon sector. 

The paper is organized as follows. In the next section, a motivation and a detailed 
description of the model is presented. The gap equation describing quark ground- 
state properties is derived. In section 3, we first describe the particular ansatz 
which produces a remarkable confinement property. The absence of quark-anti- 
quark thresholds is explicitly demonstrated for the scalar correlation function. We 
then show that this particular ansatz is indeed a solution of the gap equation. The 
phase structure of the model is discussed in some detail. In section 4, we describe the 
implication of the model on chiral properties. There is no surprise here as the model 
is designed to reproduce the chiral structure of QCD (section 2). We reproduce 
the corresponding low-energy theorems by first solving the Bethe-Salpeter equation 
for the pion field, normalize the amplitude by the electromagnetic form factor and 
extract the pion decay constant from its definition. We then establish that the Gell- 
Mann-Oakes-Renner (GMOR) relation is valid in our model. Finally, we compute 
the pseudoscalar correlation function. The sections 5 and 6 are devoted to the 
temperature and density dependence of the vacuum properties of the model. We will 



1 For SU(2) Yang-Mills theory, the chiral symmetry group is much larger, i.e. SU(2Nf), since 
the SU(2) gauge group is pseudoreal [^|. More on this later. 
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show that a deconfinement phase transition occurs at some large temperature and/or 
density. We will find that the deconfinement phase transition is accompanied by the 
restoration of the spontaneously broken chiral symmetry. The model predicts the 
two transitions at the same critical point. The final section contains some concluding 
remarks. 



2 The Model 
2.1 Motivation 



It is known for a long time that the lattice formulation of Yang-Mills theory pro- 
vides information on the confinement of quarks flnj. Subsequently quark liberation 



due to temperature was studied with the lattice version of QCD. It turned out that 
quantities which are dominated by the properties of the gluonic sector, e.g. the gluon 



condensate, vary smoothly throughout the deconfinement phase transition |25|, |26 [. 
This suggests that quark liberation may be described solely by a dynamical effect 
of the quark sector. If so, the description of this phase transition should be feasible 
by an effective quark model. We hope to gain qualitative insight into the nature 
of quark liberation by including temperature effects in the quark loops, while the 
temperature dependence of the quark interaction, mediated by the gluonic sector, 
is kept temperature-independent. 

The properties of the light pseudoscalar mesons, in particular those of the scalar 
meson and the pion, are protected by chiral symmetry of QCD and might therefore 
serve as a convenient test ground to explore other features of low energy QCD. We 
have learned in the past years that at low energies, the pion physics is phenomeno- 
logically well described in terms of an effective quark model, in which the quarks 



interact via a local current-current interaction of the NJL type |19| . However, 
the open question is: What is the signature of confinement at the level of an effective 
quark theory? 

In order to answer this question, we can be guided by a phenomenon in solid 
state physics: the localization of electrons in random potentials ||22|| , known as 
Anderson localization. In his pioneering work, Anderson showed that the key idea 
of localization can be traced back to the Hamiltonian 

H = J2(-V)a\a k + V>al 

(ik) 



where aj is the creation operator of a spin at site i of a lattice. The sum in the first 
term of ([!]) extends only over the nearest neighbors, and the constant V measures 
the strength of the nearest-neighbor hopping. This term is responsible for the spin 
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diffusion on the lattice and is the analogue of a kinetic term in continuum quantum 
field theory. The energy corresponds to a potential at site i. In the Anderson 
model, the e^'s are random variables distributed over a range —W/2 < e < W/2. A 
continuum version of this model is described by the Hamiltonian 



H = p 2 + w J2 5 ( r ~ R o 



(2) 



which describes electrons elastically scattered off impurities which are randomly 
distributed at points R a . It turns out that the precise form of the distribution of 
the potential £j is not crucial. It has been found that if the random potential exceeds 
a certain strength, i.e. 

W > W c = AKV , (3) 

where the connective constant K is characteristic of the lattice type, then the spins 
are localized at their sites, whereas for W < W c the spins are liberated , and the 
spin diffusion takes place. 

The localization of spins is intimately related to the presence of random poten- 
tials. Is there a random interaction of quarks which stems from the gluonic sector 
and which survives the low-energy limit? The answer to this question is not known. 
There is, however, a first hint from the field-strength formulation of QCD [|5|, |5|. 
When QCD (iV c = 3) is formulated in terms of the field strength, the resulting 
quark interaction takes the form 



m 



-S[T] 



cxp 



d x 



J fj, fj, 2 ^ fivJv 



(4) 



where j" = qt a ^j lJL q is the color octet current of the quarks, Tf£ = f abc T^ u is defined 
with the help of the SU(3) structure functions f abc , and V£ = {f- l )f v d p T b pv is 
the gauge potential induced by the conjugate field strength . The action of 
the field strength S[T] need not to be specified here and can be found in ^ |5|. 
Within the strong coupling limit, one observes that the gluonic vacuum decays into 
domains of constant field strength iT" giving rise to a Nambu-Jona-Lasinio type 
quark interaction in (f|). Due to gauge covariance, all orientations of the constant 
background field contribute to the quark interaction. Subsequently, the spontaneous 
breakdown of chiral symmetry was observed in the strong coupling limit 0, The 
specific form of the quark interaction due to the gluon background fields gives rise 
to a splitting of the strange- and up-quark condensates as is predicted by QCD sum 
rules |] . 

We conclude from these observations that the low energy-effective quark inter- 
action of the NJL type is in fact mediated by the gluon background field. Our key 
assumption is that in strongly fluctuating gluonic background fields (with an aver- 
age scale given by the gluon condensate), we will be endowed with a random quark 
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interaction which induces quark confinement by a mechanism similar to Anderson 
localization in solid state physics. In this paper, we shall construct a simple toy 
model which has the expected qualitative feature mentioned above of low-energy 
QCD. 

2.2 Description of the model 

In order to investigate the implications of strong random colored interactions of 
quarks, we study a model in which the quark fields are doublets of a global SU (2) 
color symmetry. We write the generating functional for mesonic Green's functions 
in Euclidean space as 

Z[<f>] = I^J VqVq expj- J d 4 x [C - q(x)(j){x)q(x)) , (5) 

G 

£ = q(x)(i(jt + im)q(x) + —[qq{x)qq{x) - qj5q(x) q^^q{x)} (6) 
+ -[qr a q{x)G aP qT l3 q{x) - q lB T a q{x) G aP q 1& T^q{x)\ , 

where m is the current quark mass. We assume that the quark interaction is given by 
a color-singlet four-fermion interaction of strength Go and a color-triplet interaction 
mediated by a positive definite matrix G af3 , which represents gluonic background 
fields. An average of all orientations O of the background field G af3 , transforming as 
G' = T GO with O being a 3 x 3 orthogonal matrix, is understood in (|5]) to restore 
global SU(2) color symmetry. Our basic assumption as motivated above is that this 
averaging amounts to assuring confinement. 

Our model is defined in Euclidean space for two reasons. First, the motivation 
of the model is provided by the Euclidean formulation of QCD. It seems reasonable 
to assume that classical Euclidean configurations such as instantons or monopoles 
might provide the random background. Second, the averaging procedure in is 
better defined in Euclidean space, since a superposition of weight functions is again a 
weight function (with a correct normalization), whereas the superposition of phases 
(as the integrand of the Minkowskian functional integral is) does not give a phase. 
The theory in Minkowski space is defined by the standard Wick rotation. We will 
have more discussions on this later. 

Instead of the current- current interaction (U), we shall use its pseudoscalar-scalar 
part that results from a Fierz transformation. The reason is as follows. As men- 
tioned, it is known that QCD in SU(2) (being pseudoreal) undergoes a symmetry 
breakdown which is quite different from what one expects in three-color QCD. To 
simulate what happens in QCD with three colors, we choose interactions so that we 
would have the correct symmetry-breaking pattern. Specifically whereas the inter- 
action of two color triplet currents which presumably follows from QCD exhibits 
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the full SU(2Nf) chiral symmetry group, the reduced interaction in @ is, however, 
only invariant under SU(Nf) x SU(Nf) transformation, as QCD is. Although we 
are dealing with an £77(2) color group, we assume that the basic idea of the confine- 
ment mechanism, developed below, does not depend qualitatively on the color group 
under investigations^]. In order to keep contact with QCD as closely as possible, we 
therefore choose our low-energy effective quark theory to exhibit the chiral patterns 
of QCD in four dimensions. 

In order to make contact with more familiar formulations of effective quark mod- 
els, we first study the limit where the color-triplet interaction G af3 is weak. In this 
limit we can perform the average over the background orientation O using a cumu- 
lant expansion. The colored part of the quark interaction becomes 



J dX /(A) exp i — — J d^x [qr a q(x) qr c 



q(x) - qi<oT a q(x)q^T a q(x)} + 0(A 2 



(7) 

where A is an eigenvalue of the matrix G af3 , and / (A) is the corresponding eigenvalue 
density. At lowest order of the cumulant expansion, we obtain the familiar Nambu- 
Jona-Lasinio model [18| with a global color symmetry. Terms multiplied by A 2 are 
eight-quark interactions. The average over the background field G al3 in (|5|) obviously 
incorporates the interaction of more than four quarks. If these interactions are not 
small or equivalently if the background fields in (|5[) are not weak, one then has 
to abandon the cumulant expansion, and instead study the quark theory of fl5|) 
with fixed background field in a certain approximation, which we will specify below, 
and average over the background fields afterwards. If the approximation used is 
not bound to weak couplings, this approach can be applied even if the cumulant 
approximation fails. 

Specifically, the approximation under investigation is to introduce meson fields 
on top of the scalar condensate of the vacuum and to treat their interactions per- 
turbatively. This approximation does not resort to small couplings and is expected 
to give good results if the number of mesonic degrees of freedom is large. This 
approximation is the usual one applied to study the physics of light hadrons in the 
context of the Nambu-Jona-Lasinio model. We will not further question the validity 
of the approximation, but investigate the ground state and the properties of the 
light mesons within this scheme. The quark interaction in @ is linearized by means 
of color-triplet (a a , n a ) and color-singlet mesons (a, tt). Integrating out the quark 
fields in the Hubbard- Stratonovich formalism, the resulting effective meson theory 



2 We do not expect the answer to the question as to whether the quarks are confined or not to 
depend on the color group used, whereas the actual value of e.g. the pion decay constant might 
depend on whether we have the SU(3) or the SU(2) gauge group. 
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IS 



C M = ~ In (i$ + iM + ia + iM a r a + ia a r a + 7r 75 + vr a r a 7 5 ) (9) 



+ 7^T ((<* + M - m - i<J) ) 2 + (tt + 5 ) 2 

where M a and M are, respectively, the color-triplet and color-singlet constituent 
quark masses, and we have decomposed the external source into a scalar and 
pseudoscalar parts, i.e. = 0o + 057s- The approximation mentioned above consists 
of truncating the expansion of (^) in terms of the meson fields. The zeroth-order 
approximation provides access to ground state properties, in particular, the quark 
condensate. From 

5\nZ[0] 5lnZ[0] 



0, . 1 J = (10) 



Sa a ' So- 



one obtains 



- — Tr I . 1 + — (M - m) = , (11) 

--Tr{ ' B A + fa- 1 )* tin = , (12) 

where V4 is the Euclidean space-time volume, and the trace extends over internal 
degrees of freedom as well as over space-time. In the latter case, a regularization 
is required, which we will specify when needed. Different solutions of the equations 
(O-T^) correspond to different phases of the system. 



3 The Confining Phase 



Here we will show that a particular solution of the gap equations (JTT|) — (^2|) with 
an imaginary color-triplet constituent mass exists, i.e. M a = —iM a with M a real. 
Before discussing in detail the existence of such a solution, we first illustrate its 
remarkable physical consequence. 



S 



3.1 The scalar correlation function 



Consider the color-singlet scalar correlation function 

A s (p) = J d 4 x e-**" (qq(x) qq(0)) , (13) 
within the effective meson theory (P). Its connected part is given by 

J d0o(x)d0 o (O) 

To compute this, one has to consider the fluctuations of the scalar color-singlet field 
o"o as well as those of the fields a a , since, for a fixed interaction G Q/3 , the colored 
mesons couple to er by a quark loop (see the first term on the right-hand side of 
(0)). Fluctuations of the pion fields n a , n do not contribute to the scalar correlation 
function, since the pion fields have the wrong quantum numbers. Expanding (P) up 
to second order in the meson fields, the action, in momentum space representation, 
is 

S (2) = {l^(p)Tl s (p>(-p) + l^(p)Uf(p^(~p) (15) 



where 



+ ta a (p)K a (p 2 )a(-p) - _L0 o (p)0 o (- p ) - J-^ai-p) 



IW) = - J -j^tr {S{k + p)S{k)} , (16) 



r d k 

nf (P 2 ) = (G-'r" - J j^tr{r a S(k+pyS(k)} , (17) 
r d A k 

K a (p 2 ) = -i ——ti{r a S(k + p)S{k)}. (18) 



„ (2tt)* 

The quark propagator S(k), in momentum space, is 

S ^ := $ + t(M + tM«r«) ' (19) 

Our model @ is defined with an average over all global 577(2) orientations of the 
interaction matrix G a @. In order to render the averaging procedure easy, we integrate 
out the colored mesons in (|I5D , i.e. 

s$, = / (0i {r(p) [ n ° + ^(n; 1 )^] *(- P ) (20) 

1 i 

-Mp)M~p) - TrMp)°{-p) 
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(2) 

In fact, we will observe that S e ^ no longer depends on the orientation of the inter- 
action matrix, so the averaging procedure is trivial. 

In order to study the effective theory of the singlet meson a, we explicitly calcu- 
late the polarization functions 11°, 11"^ in flHf ) and (|H]), respectively, as well as the 
mixing K a flUf). For this purpose, it is convenient to introduce the eigenvectors of 
the matrix M a r a , i.e. 

M a r a \±) = ±M\±), M v.\ I . (21) 

They possess the property 

M a 

(±\r a \±) = ±— , (22) 

which will be extensively used later. The detail of the calculation is left to Appendix 
[A|. It turns out that the quantities ([16]) and ([18]) can be expressed in terms of two 
functions H Q (p 2 ), H v (p 2 ), i.e. 

1 M a 
IW) = £T " MP 2 ) , #V) = — H v (p 2 ) . (23) 

For later convenience, we explicitly write them down here: 

h ^ - 4 W!H^ + < m -- m >}- (24) 

H.tf) = -U £ ia J ^ - (M -M)} (25) 

where 

Q = a(l- a)p 2 + (M + iM) 2 . (26) 
We also find (see Appendix |A]) that M a is an eigenvector of the polarization matrix 

s M \G C J M K J 

where we have used the property (of which we will have more to say in the next 
subsection) of the solution of the Dyson-Schwinger equations ( ]IT| , [T2"| ), namely, that 
M a is an eigenvector of the symmetric matrix G a/3 with eigenvalue G c . It is now 

(q\ _____ 

straightforward to calculate S^/f, (EO). Since K a is proportional to M a and using 



that the eigenvectors of the symmetric matrix G a ° are orthogonal, one obtains 



G a * ' ±-H lt (tf) 



a(-p) (28) 



1 i 

-Mp)M~p) - ttMp)°{-p) 
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Note that SfL depends only on M and M, which are invariant under SU{2) rota- 
tions. This is the desired result, since the average over the SU(2) orientations can 
now be trivially performed. 

We are now going to study the occurrence of an imaginary part of the scalar 
correlation function signaling a quark- ant i- quark threshold. The crucial observation 
will be that whenever the trace of the color indices is required, the contributions 
with M a occur in conjugate complex pairs. In particular, this will erase imaginary 
parts of the correlations function, and no quark-anti-quark threshold will occur. 

In order to work out this phenomena in some detail for the scalar correlation 
function (|T3|) , it is sufficient to study the functions H (p 2 ), (|2~4f) and H v (p 2 ), fl25[) . 
since they provide the complete correlator with the help of (p8[). We rewrite, for 
instance, H (p 2 ) as 

= ht^C^ f 1 - ZTq + ¥Tw] + (M - - M) ' (29) 

To illustrate the disappearance of the quark-anti-quark threshold for M 7^ 0, we 
first study its occurrence for M = 0. In this case our model describes the scalar 
correlation function in the usual constituent quark model with a constituent quark 
mass M . The term of interest is the second one in the curly bracket in fl2Tj|). After 
integration over u, this term becomes essentially InQ. This implies that whenever 
Q becomes negative, the function H , (|2~9"D , acquires an imaginary part. In order for 
Q to become negative, the Euclidean momentum p 2 must satisfy, 

-p 2 < 4M 2 , (30) 

implying that the quark-anti-quark threshold occurs at a (Minkowskian) momentum 
Pm = 2Mo, which is the familiar result. For M = the function H v (p 2 ) does no 
harm, since it is identically zero by definition, (E^). 

We are now going to show that for M 7^ no threshold will occur at all. Adding 
up the contributions from M and — M in Q2*9|), the crucial term becomes 



™i 3 f 1 j f j (u + W)W +4M 2 M 2 

where W = a(l — a)p 2 + Mq — M 2 . An analogous result holds for the function 
H v (p 2 ). One finds 

3 r 1 r A2 ?M 2 M 2 

m rit = ~ da du -. ^ 9 - . (32) 

2tt 2 Jo Jo [u + W} 2 + 4M 2 M 2 V ; 

We find that the logarithmic divergence is screened if MM 7^ 0. No imaginary part 
occurs for M 7^ and M/0. This is our main observation. For non-vanishing 
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current mass m, one always expects at least a small constituent quark mass M . 
Therefore, the main ingredient in avoiding the decay of the scalar meson into a 
quark-anti-quark pair is the non-vanishing value of M. In the following, we refer to 
the phase of the constituent quark model (B|) with M / as confining phase. 

In the chiral limit, the chiral symmetric phase (Mo = 0) needs further discus- 
sion. In section (|^) we will find that temperature induces a deconfinement phase 
transition, and that chiral symmetry is restored at the same time. Here we find 
by an inspection of (0,^) that the restoration of chiral symmetry (M = 0) is 
accompanied by the occurrence of quark-anti-quark thresholds. 

3.2 Phase structure 



Here we will search for solutions of the gap equations ( jiT| - |T2|) with a non- 
vanishing, imaginary constituent quark mass in the color-triplet channel, which im- 
plies the remarkable consequences discussed above. We will discuss the dependence 
of the phase structure, and, in particular, the phase transition from the confining 
phase (M ^ 0) to a non-confining phase (M = 0), on the parameters of the model, 
i.e. Go, G a ^ and m. For this purpose, we have to analyze the solution of the gap 
equations (p"T|-|T2|). In order to solve fll2|), we assume M a to be an eigenvector of the 
matrix G af3 , i.e. 

G^M 13 = G c M a . (33) 
This reduces eqs. ([□]) and (|T2"D to (see Appendix ||) 

-L(Mo-m) = -M I (M ,M) + MI V (M ,M) , (34) 



G 



o 



i-Af = -MI (M ,M) - M I V (M ,M) , (35) 



where 



1 r A2 u + M 2 — M 2 

/o(M °' M) = -2W0 dU U + 2u(M$ - J 2 ) + (Mq + ' (36) 

1 f A2 M M 
4(Mo ' M) = T*h dUU u* + MM - M 2 ) + (Ml + M 2 ) 2 ' (37) 

where a sharp 0(4) cutoff A was introduced. Since I v is proportional to M M, 
M = is always a solution of eq. (^). In this case, eq. (51) is reduced to the 



gap equation of the standard Nambu-Jona-Lasinio model with an additional SU(2) 
degree of freedom. This implies that the standard Nambu-Jona-Lasinio model is 
contained in the extended model (|J) as a special case, and we expect the theory (||) 
to be phenomenologically as successful as the Nambu-Jona-Lasinio model. 
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The solution of the system of eqs. (|3^j35|) was investigated numerically. Figure 
[I] shows the color-singlet and color-triplet constituent quark masses, M and M, as 
function of the color-singlet coupling Gq. For sufficiently small color-singlet coupling 
strength a confining phase with M^O exists. Also shown is the non-confining phase 
with M = for all values of Gq. In order to decide which solution forms the vacuum, 
one has to compare the classical action, i.e. 

Ac = J (0 {-2 H(k> - Ml - M 2 ) 2 + 4M 2 A; 2 ] - ^M 2 + ±-(M - m) 2 

(38) 

of both solutions. One finds that the confining solution has a lower action and forms 
the vacuum. 

It is observed that the color-triplet coupling strength G c must exceed a critical 
value in order to allow for the desired imaginary constituent quark mass in the 
color-triplet channel. Figure |2| shows the critical coupling as function of the color- 
singlet coupling Go for different values of the current mass m. A strong color-singlet 
constituent mass Mq seems to suppress the occurrence of the confining phase. The 
situation can be compared with that of free electrons in a solid. It was discovered 
by Anderson that the electrons get localized if the density of impurities exceeds a 
critical limit |2^| . 

Finally, we present the result for the scalar correlation function A c s (p 2 ), ( |i4|) 
which we have calculated in the last subsection. The final result (two integrations 
are left to a numerical calculation) is given by 

A c (p ) = -pr — 7^2 ~j . 9 , ; h?(v 2 ) ' 
Go G J- _ H (p 2 ) + 

Figure |] shows the correlation function A c s {p 2 ) as function of the Euclidean momen- 
tum transfer p 2 for different values for the current quark mass m. From Figure 
we conclude that increasing the current mass m at fixed values of Go and G c drives 
the system towards the deconfinement phase transition. For large m (close to the 
critical value of m where the transition to the deconfined phase occurs) one observes 
a resonance-like peak at a negative momentum squared, which is reminiscent of the 
quark anti-quark threshold of the deconfined phase (M = 0). In the latter case, 
our model is identical to the Nambu-Jona-Lasinio model. It is well known that 
in the standard Nambu-Jona-Lasinio model, a weakly bound scalar particle occurs 
in the spectrum. It has a mass of two times the constituent quark mass implying 
that the particle pole occurs at the threshold position. The physics of the scalar 
meson is therefore beyond the scope of the standard Nambu-Jona-Lasinio model. 
This defect of the NJL- model is cured in our model (0): the threshold is avoided 
in the confining phase (M ^ 0). At the same time, no scalar particle is present 
in the spectrum. There appears a "resonance" near the threshold position with a 
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width proportional to M. This result is in agreement with nature, since no scalar 
bound state is observed in the meson spectrum at the relevant energy scale. The 
scalar Green's function develops a structure around the threshold position without 
allowing the meson to decay into two quarks. This might be a precursor to a parton 
structure of hadrons at high energy. Unfortunately we cannot push this conjecture 
further, since our model is limited to low energies. 

One might raise the question as to whether this result is in not disagreement 
with dispersion relations, since the correlation function we have is a real function 
in the entire momentum space. It is interesting to note that our model does not 
allow a scalar particle to appear, although a scalar field is needed in the bosonization 
approach to represent the quark interaction. The scalar correlation function at hand 
is, therefore, not appropriate for discussing dispersion relations. On the other hand, 
as we will see soon, the pion-pion correlation function exhibits the pion pole and is 
otherwise real, so one might worry about the dispersion relation in this case. Note, 
however, that in order to arrive at this result, we have expanded the effective meson 
theory (Q) up to second order in the meson fields ignoring mesonic interactions. If 
these interactions are included, we expect to recover the physical thresholds of an 
interacting meson theory satisfying dispersion relations. 

Although the correlation function is real, a non-trivial structure appears at the 
would-be threshold position. A close inspection of the correlator (|39| ) (see Appendix 
|CD shows that the function possesses cuts in the complex p-plane. One might then 
question whether these complex structures violate fundamental principles of quan- 
tum field theory (e.g. stability). We note that our approach - the bosonization 
procedure of section and the related approximations - is supposed not to violate 
any of the analyticity requirements. In comparison to other NJL-type models, our 
ground state (vacuum) is simply more complicated, with the standard NJL model 
corresponding to a particular case (M = 0). It may be that some of the fundamental 
requirements of QCD are not correctly implemented. If so, then our detailed calcu- 
lation may be providing mechanisms for the cancellation of unwanted features. In 
this case, we could expect our model to give some insight in the analyticity structure 
of Green's functions of a confining theory. This question is of particular interest, but 
is beyond the scope of the present paper. See for a discussion on this matter. 
To check that nothing is amiss with our model, we have studied the analytic struc- 
ture of the scalar correlation function in Appendix |C| and explicitly verified that the 
stability criterion is indeed satisfied thanks to a cancellation mechanism mentioned 
above. 
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4 Chiral Properties 



In this section we study the properties of the pion that should emerge as a 
Goldstone boson of the spontaneously broken chiral symmetry. We will show that 
our model is in agreement with the low-energy theorems, and will verify by an 
explicit calculation that the Gell-Mann-Oakes-Renner relation holds. 



4.1 The pion Bethe-Salpeter equation 



One way to extract the properties of the pion is to study the pseudoscalar cor- 
relation function in the color-singlet channel which can be obtained in the same 
way as for the scalar correlator discussed in section |3.1j . Here we prefer to use a 
different approach which is to calculate the pion Bethe-Salpeter amplitude. This 
method will illustrate the role of the hidden color components of the pion, and will 
readily provide an access to such observables as the pion decay constant and the 
pion electromagnetic form factor. 

The Bethe-Salpeter equation for the pion amplitude (P , P a ) is directly obtained 
from the effective meson Lagrangian Cm (§) by 

For a fixed orientation of the interaction matrix G a ^ , the left-hand side of this 
equation becomes 

/ £ + Tr{ l5 S(k + ph 5 S(k)} Tr{T« l5 S(k + p) l5 S(k)} W P \ 

{ Tr{r a l5 S(k + p) l5 S(k)} (G' 1 )^ + Tr{r a -f 5 S(k+p)r l3 -f 5 S(k)} J { P a J ' 

(41) 

where the trace Tr extends over the momentum space (k) as well as Lorentz- and 
color-space. One observes that the ansatz 

M a 

P a = i^Pi (42) 
for the color-components of the pionQ reduces ([II]) to (see Appendix 0) 

{ Up") ^+Up 2 ) J [ Pi(p 2 ) J " u ' {A6) 



Once the average over all orientations of the int erac tion matrix is performed, these parts 



become the hidden-color components (compare section 4.3). 
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where the functions X and X v are defined by 

<r / ?s 1 f 1 r f A \ u - a(l - a)p 2 + A 2 /w 

j Qo 2 = -/ rfa / duu i - n + (M^-M), (44 

; 4tt 2 Jo To u + a (l - a)p 2 + A 2 ] 2 V ; ' V 1 



-r > i f 1 , f A2 , u - a(l - a)p 2 + A 2 

UP 2 ) = j- 2 da / rfwn-— -A ^——5 - M^-M , (45) 

47T Jo Jo [u + a(l - a)p 2 + A 2 ] 



with A = Mo + «M. An inspection of (|4q) yields that both Xd and X„ are real. 
Demanding (|43| ) to have a non-trivial solution for P 2 = — m 2 leads to a nonlinear 
equation to determine m n . Instead of solving this equation numerically, it is more 
instructive to solve it analytically for pion masses small compared with the con- 
stituent quark mass M . For this purpose we expand the functions X and X„ up to 
linear order in p 2 , i.e. 

X (p 2 ) = Jo + F p 2 +0(p 4 ) , l v (p 2 ) = I v + F v p 2 + 0(p 4 ) , (46) 

where J and I v depend, respectively, on M and M (see fl5E|) and |57D), and the 
quantities F , F v are defined by 

otv z Jo [k 2 + (M + «M) 2 ] 

F " = "A r ^) k " r?2 , 7^- .^ 21 2 -IM-+ -M) , (48) 
87r 2 Jo [k 2 + (M + iM) 2 } 

A non-trivial solution for the Bethe-Salpeter amplitudes Po, Pi exists provided 

With the help of the equations of motion (|34l) and (|35|) , i.e. 

1 M m 1 M 

Go~ + / ° " M^ + MoGo"' _ ~~M ' (50) 

eq. ( fi^D can be rewritten as (setting p 2 = —m 2 ,) 

mlf 2 c = Am (qq) + O(mJ) + 0(m 2 ) , (51) 

where we have used that the quark condensate (qq) is given by M /G , and where 
f 2 c is defined by 

f 2 c = 4{m 2 F - M 2 F - 2MM F,} . (52) 

Equation fl5"T|) is the Gell-Mann-Oakes-Renner relation. It tells us that in the chiral 
limt (m = 0) the pion is massless if chiral symmetry is spontaneously broken ((qq) ^ 
0). Below, we will show that f 2 c , defined in ([52]), coincides with the pion decay 
constant in the chiral limit (m = 0). 
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4.2 The electromagnetic form factor 



Since the Bethe-Salpeter equation provides only the relative weights of the am- 
plitudes Pq, P a , an additional physical input is needed for the normalization of the 
amplitudes. One may use the electromagnetic form factor of the pion. The elec- 
tromagnetic form factor provides the desired additional constraint. The form factor 
F(q 2 ) is defined via the electromagnetic vertex [J27|, i.e. 



<7r(f/)|J^(0)|7r(p)> = ie(p'^+p,)F(q 2 ) } (53) 

where is the electromagnetic current and q = p' — p. In order to normalize the 
pion charge to unity, we demand that 

F(0) = 1 . (54) 

Once the form factor is calculated, the (mean-square) charge radius of the pion is 
simply given by 

2 / 2\ R 9F(q 2 ) 

r ms = ( r h = - 6 o 9 U 2 =o > ( 55 ) 



where q is the momentum in Minkowski space. Writing the electromagnetic current 
in terms of the quark fields, 

J fi( x ) = Q{ x )ln<l{ x ) > ( 56 ) 
we can study the matrix element fl5"3"|) most efficiently by using the Bethe-Salpeter 



amplitude of the pion. Its graphical representation is given in Fig.|Pj. For p' — 0, 
the matrix element is 

d A k f / M a aT ^\ an p. an p. 



/^tr{ 75 ( Po + i ^r« Pl )s { k- P - M k + l) (57) 

-r^Pi) S(k - f ) } . 



/ M 13 
15 \ M 



The evaluation of this matrix element is left to Appendix [E]. With (|E.7| ), the form 
factor, for small momentum, reads 

F(p 2 ) = (Pg-P?) (F -p 2 R ) - 2P Pi (F v -p 2 R v ) . (58) 



4 Since the diagrams (in particular, (a)) are divergent, the result depends on the actual choice 
of the loop momentum. A more sophisticated cutoff procedure (e.g. Schwinger's proper-time 
regularization) would remove this ambiguity. 
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From this expression follow the desired normalization ( [54]) of the pion Bethe-Salpeter 
amplitude and the pion charge radius (p3|): 



We can now discuss the full momentum dependence of the form factor F(p 2 ). 
The expression for F(p 2 ) is derived in Appendix |]. The integration over the angle 
variable and the radial component of the momentum space was performed numeri- 
cally. The numerical result is presented in Figure [5] There exists no vector-meson 
poles at negative momentum transfer, but instead a peak, its width strongly de- 
pending on the parameters of the model. Presumably the resonance-like structure 
is an artifact of the model, unrelated to any physical processes. It probably reflects 
the suppression of the quark-anti-quark pole present in the standard NJL model 
without confinement. As mentioned, the physics of the vector mesons is beyond the 
scope of a local four-quark interaction^]. We expect that a more realistic model with 
quantum loops will produce a vector- meson (i.e., p) pole instead of the resonance. In 
fact, we observe a pion pole and no further structure in the pseudoscalar correlation 
function in our model, which is supposed to give good results for the pion physics 
(see below). 

It is quite remarkable that we obtain a real form factor for arbitrary large 
(Minkowskian) momentum manifesting the absence of quark and anti-quark thresh- 
olds. This cures an outstanding problem of the standard non-confining NJL-model. 
This may have an important consequence in heavy-meson physics. For instance, 
the non-perturbative description of the electroweak currents of a heavy and a light 
quark is beyond the reach of the usual NJL-model due to the occurrence of quark 
anti-quark threshold effects. Electroweak currents such as the one considered here 
enter into the decay rate of a B-meson into a pion and an electron. In the limit of 
the electron being very fast, one should be able to extract the element Vb u of the 
Kobayashi-Maskawa mixing matrix from the differential cross section. However, it 
has been shown that in the limit of a fast electron, non-perturbative contributions 
to the electroweak currents become important. We believe that a more realistic 
formulation (SU(3) color) of our model might be able to provide the desired results. 

The normalization fl59|) completes the calculation of the pion Bethe-Salpeter 
amplitude. This is the desired result, since we are now able to calculate matrix 
elments involving pions. In particular, we will obtain the pion decay constant in the 
next chapter by a direct calculation. 



5 In particular, we do not have any vector and axial-vector channels in the quark interaction 
, which are known to be important for the physics of vector mesons. 




(59) 
(60) 
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4.3 The pion decay constant 



An expression of pion decay constant was already deduced in section |4.1| by de- 
manding that the Gell-Mann-Renner-Oakes relation ([51]) be quantitatively satisfied. 
The aim of this subsection is to show that the decay constant (0), obtained there, 
is identical to that derived from its definition. 

The pion decay constant /„. is defined by the matrix element 



(0\A,(x)\n(p))\ p ^ = iJv* V »U 



(61) 



which describes the coupling of the pion to the vacuum via the axial vector current 



Afi,(x) = ( l{x) r Y5 r Yfj,q{x). The matrix element (|B~1]) , also shown in Figure |], is 



(62) 



An average over all color orientations of the interaction matrix G af3 - which is 
equivalent to an average over all color directions of M a - is understood in (|62|). A 
straightforward calculation of for a given vector M a yields 



d A k 



(27T) 



.4 



[k 2 + A 2 ]' 



+ (M — > —M) 



Pa 



(63) 



H) 



d A k 



A 



(2vr) 4 [[k 2 + A 2 Y 



(M -> —M) 



Pi 



This result can be further simplified by introducing the functions Fq and F v from 
(|^,||), i.e. 

2i Vil {(M F - MF V ) P - (M F V + MF ) P x ) . (64) 
The ratio e of P over -Pi is provided by the Bethe-Salpeter equation (|43"D, i.e. 



Pi 
Po 



G 



l v (p 2 = -ml) 
T + J o(p 2 = -ml] 



(65) 



whereas the overall normalization is constrained by the electromagnetic form factor 
(|59|). Expressing P\ in terms of e and Po and eliminating Po with the constraint 
(m), we obtain the final result 



J TT 



[M F -MF V - (M F v + MF )e] 
(l-e 2 )Fn- 2eF, 



(66) 



where the functions Po, F v are given by (0) and (|4"5D. This is our result for the 
pion decay constant, valid for all values of the current mass. This decay constant 
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needs to agree with the one extracted from the Gell-Mann-Renner-Oakes relation, 
(0), only in the chiral limit (m = 0). Before we show that this is indeed the case, 
we would like to comment on the result fl66|). First note that the average of all 
orientations has still to be performed. Since our result (|5oT), however, depends only 
on the invariant quantity M rather than on M a , this average is trivial. Note further 
that a term proportional to Pi enters into the physical decay constant. This is the 
contribution from the hidden-color components of the pion. 

In order to calculate /„- in the chiral limit, we proceed with (|63|). From the 



Bethe-Salpeter equation (f23|) in the chiral limit p = —m^ = 0, one finds 

With this result, the decay constant in the chiral limit becomes 

U c = 2{M 2 F - 2MM F V - M 2 F } g . (68) 

Using ( |S7D in the normalization of the Bethe-Salpeter amplitude (|59"D , one may 
eliminate P\ in (|68|) to obtain 

fl c = 4{(M 2 -M 2 )F - 2MM F v j . (69) 



This expression is identical to (p2|). This establishes that the Gell-Mann-Oakes- 
Renner relation is indeed valid in our model. 

We finally present the numerical data for the pion mass and the pion decay 
constant. We have numerically solved the Bethe-Salpeter equation (f43|) to obtain 
the mass of the pion and the ratio e = Pi/Pq. The decay constant /„. was then 
calculated from fl6"6"|). The result for m n and as a function of the current mass 
is shown in Figure One observes that m 2 depends almost linearly on the current 
mass with the slope given by the Gell-Mann-Oakes-Renner relation. The decay 
constant decreases with increasing current mass. 

4.4 The pseudoscalar correlation function 



It is instructive to compare the result of the scalar correlation function with that 
for the pseudoscalar correlator where one expects a pion pole. The pseudoscalar 
correlation function is defined by 

J d(j) 5 {x)d(j) 5 {0) 

In order to derive A£, we expand the effective meson theory @ up to second order 
in pion fields. Since an average over all orientations of the interaction matrix G a @ is 
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required, it is convenient to integrate out the colored pion fields. One observes that 
the resulting theory for the color-singlet pion does not depend on this orientation, 
so the averaging is trivial. Since the explicit calculation parallels closely that of the 
scalar correlator in section pO, we shall simply present the final result: 



S { ;/ f = J {-rr(p)Tl 5 (p 2 )n(-p) + _0 5 ( p )0 5 (- p ) + —fc^-p) 

where the pion dispersion formula U 5 (p 2 ) is given by 



(71) 



IW) = + W) + i Ir 2 , ■ (72) 



The functions Xo and X v were defined in (|44] ) and ( |4q) , respectively. The dispersion 
formula II5 is shown in Figure |7| as function of the Euclidean momentum transfer in 
the chiral limit (m — 0). Note that the function ^(p 2 ) is zero at p 2 = 0. This zero 
gives rise to a pole of the correlation function at zero momentum transfer confirming 
the pion as Goldstone boson. The most striking feature of Figure [7| is a singularity 
at positive Euclidean momentum. It arises from the second term in ([72]), which 
stems from the integration over the colored pions. The singularity occurs at the 
momentum where the colored pions go on-shell. The occurrence of the singularity 
therefore is quite natural and does not depend on the details of the model. In 
our model, the colored pions go on-shell at positive Euclidean momentum p 2 , which 
might indicate that the colored pions condense. In this case our approach, expanding 
the Lagrangian up to second order in the colored pion fields and integrating them 
out, is no longer appropriate. We will leave this issue to future investigations and 
here study the pion pole in some detail. 

One can deduce the pion mass m 2 from the position of the pion pole in the 
correlation function A^(p 2 ). It is worthwhile to check whether this result for the 
pion mass agrees with that obtained from solving the Bethe-Salpeter equation (see 
section 4.2). One observes that both masses are indeed identical, because the con- 



straint for the dispersion relation to be satisfied coincides with the condition ( |49|) 
which guarantees a non-zero Bethe-Salpeter amplitude. Finally the numerical result 
for the correlator (|70| ) is presented in Figure ^. The correlation function is entirely 



dominated by the pion pole. This is compared with the result for the scalar corre- 
lator. In the latter case, no pole occurs at all, but the scalar correlation function is 
significantly influenced by a peak structure. 

A second pole in the correlation function occurs at a positive value of the Eu- 
clidean momentum p 2 . This pole is due to the fact that the occurrence of the colored 
pion pole in the dispersion relation il 5 gives rise to a further zero of U 5 (p 2 ). The 
pole of A£ at a positive momentum squared is intimately related to the properties 
of the hidden-color pions and will be the subject of a future work. 
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5 Temperature Effects 



The model (||) offers a confining phase and a deconfining phase, both emerg- 
ing from the classical equations of motion (ITT 



In the section given above, we 
investigated the phase structure in the model's parameter space. We found that 
the ground state is in the confining phase if the color-triplet coupling strength is 
strong enough. In this section, we will investigate the influence of temperature on 
the phase structure. In particular, we will start from the system with the vacuum 
in the confining phase and will study the type of phase transitions, if any, to the 
deconfining phase. 

For this purpose, we must generalize the gap equations ( pl - |i~2"| ) to finite temper- 
ature. The trace in the first terms on the left-hand side stems from the integra- 
tion over the quark fields. In order to introduce temperature, we adopt the usual 



imaginary-time formalism [2"4~] and confine the configuration space of the fermionic 
fields to the configurations which are anti-periodic in the Euclidean time direction 
with a periodic length 1/T with T the temperature. At finite temperature, the 
integration over the zeroth component of the Euclidean momentum in the terms of 



( |i~I}]T2| ) is replaced by a discrete sum over Matsubara frequences. In order to illus- 
trate the evaluation of such trace terms at finite temperature, we explicitly work 
out a term, which arises in (|TT]-[T2|) once the color trace is performed, i.e. 



F := Tr 



k 1 + i(M + iM) J 

Performing the Dirac trace as well as the trace over space-time, one obtains 

d 3 k M + iM 



(73) 



(2tt) 3 ix 2 T 2 (2n + l) 2 + (a + ib) 2 
where V is the space volume, and the quantities 



(74) 




- k 2 + M 2 — M 2 + y(k 2 + M 2 - M 2 ) 2 + 4M 2 M 2 



(75) 
(76) 



were introduced as an abbreviation. In order to split the zero temperature part, 
which is cutoff-dependent, from the temperature-dependent contributions, which 
are finite, we perform a Poisson resummation of the Matsubara sum in ([74]), i.e. 



AV / dn 



d 3 k 



M + iM 



(2tt) 3 n 2 T 2 (2n + l) 2 + k 2 + (M + iMf 



F{y ± 0) . (77) 
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The first term is the term with zero conjugate Matsubara frequency [v — 0), which 
yields the zero temperature result if we substitute fc — ^T(2n + 1) for n in the 
integral in (|77|) . It is this term which needs regularization. It will be performed by 
introducing a sharp 0(4) cutoff as it was done in the zero-temperature case. The part 
of the non-zero conjugate Matsubara indices is finite and needs no regularization. 

F( ^ 0) ^^Mi±lE [e ^ + l] -\ (78) 

1 a + ib 1 J 
Evaluating the trace terms in ( |Il"| - |T2|) as sketched above, one finds that the gap 



equations acquire an additional term that contains a temperature dependence. In- 
troducing 

ae a ' T cos(6/T) + a - be a l T sin(6/T) 
R ~ {a 2 + b 2 ){e 2a / T + 2 cos(6/T)e a /* + 1) ' ( ' 

ae a ' T sin(6/T) + b + be a ' T cos(6/T) 
1 ~ {a 2 + b 2 ){e 2a / T + 2 cos(6/T)e a /* + 1) ' ( ' 

the modified gap equations are 

G a \M -m) = M I + (M ,M) - 4r I dk k 2 (M S R + ¥S;) . (81) 
G~ X M = M/_(M ,M) - 4t fdkk 2 (MS fi - M S 7 ) , (82) 

TV 2 J 

where the functions I± are defined with the help of (0,[37]) by 

7 + = + wj- • l - = - TT 1 " • (83) 

This is the main result of this section; the temperature-dependent color-singlet and 
color-triplet constituent quark masses will emerge from these equations. 

We have studied numerically the solutions M and Mo of (^T]-^) as function of 
temperature. The coupling strengths Go and G were chosen in order for the system 
to be in the confining phase at zero temperature. The result is shown in Figure 
([|). If the temperature exceeds a critical value, a first order phase transition from 
the confining phase (M ^ 0) to the deconfined phase (M = 0) takes place. The 
deconfining phase transition is accompanied by a sudden drop of the color-singlet 
constituent mass M indicating at the same time the restoration of chiral symmetry. 
The small residual constituent quark mass is due to the current mass m = 0.02A 
which explicitly breaks chiral symmetry. 

For phenomenological applications, the dependence of the deconfining phase 



transition on the current quark mass is of particular interest. Figure 10 shows the 



critical color-triplet coupling strength G c as function of the temperature for some 
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values of the current mass. The result suggests that quark liberation of all quark 
flavors occurs approximately at the same temperature (assuming that the system 
was in the confining phase at zero temperature). 



6 Finite-Density Effects 



In the last section, we observed a phase transition from the confined phase to the 
deconfined phase at high temperature. For phenomenological reasons f28|, one also 
expects a phase transition to occur at high baryonic density. Here we investigate 
the existence of this phase transition in the random background quark model fl5|). 

In order to study the system at a non-zero baryonic density, we introduce a 
chemical potential fi into the Lagrangian (||): 

£d = q{x)(i$ + im + iftj )q(x) + G [qq(x) qq(x) - qj 5 q(x) qj 5 q(x)) (84) 
+ [qr a q(x) G a ^qr p q(x) - q^T a q{x) G af3 q^q{x)} , 

The bosonization procedure - introducing scalar and pion fields - is unchanged 
by the presence of the chemical potential. The gap equations (|TT1, 12) acquire an 
additional part due to finite density, i.e. 





Tr 






where 




h 




I 13 


= l T r 

v 4 



icf) + iM Q + iM a r 
i 



I F (M ,M) + i-(M Q 



m) 



i$ + iM + iM a r c 



-T 



P 



4(M ,M) + (G^Y'Mi 



0, (85) 
), (86) 



i$ + ifi'jo + iM + iM a r a 
i 

i$ + i>7o + iM + iM a r a 



V, 



-Tr 



l T r 

V 4 



i@ + iM + iM a r a J 
' i 
ity + iM + iM a r c 



"0 



One easily verifies that both functions Ip and I F are finite and need no regular- 
ization. In order to calculate these functions, we first perform the color trace by 
introducing the eigenvectors of the matrix M a r a defined in (PH). After taking the 
trace over Dirac indices, the function Ip becomes 



I F {M ,M) 



d 4 k 



4(M + iM) 
(fco + ^) 2 + (« + *) 2 ~ H + (a + ib) 2 



4(M + iM) 



+ (M -> —M) 



i9) 
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where the functions a and b depend on the momentum k as defined in (|75|j76| ). It is 
now straightforward to evaluate the k integration in (|55|): 

I F (M , M) = -i f M + (M - -M) . (90) 

Ja(k 2 )<n (2-k) 6 la-b 

One observes that the color trace again enforces that the function Ip be real. Anal- 
ogous considerations hold for the function I F . The final result for the gap equations 

G \M -m) = M I + (M ,M) - \ dk k 2 M °° + f & , (91) 

tx a Jo cr + cr 

G-H4 = M/_(M„,M) - Lj k ' dk e'^f^, (92) 



where the functions I± are defined in (83) and a(k 2 ) and b(k 2 ) are defined in ([75]) 
and (f76|) respectively. The upper bound of the integration over the three-momentum 
k is provided by the Fermi sphere of radius kf where kf is defined by a(k 2 ) = //. 
The solutions of the coupled system (OT|J2|) provide the constituent quark mass M 



and the mass M in the colored channel as a function of the chemical potential \x. 
For physical applications, it is convenient to express the chemical potential in terms 
of the baryonic density defined by 

(d\TiZ , . dlnZ, A , x 



From (|5|) and the Lagrangian Lp, at finite chemical potential fi in (S3), one obtains 

d\nZ / f -1 1 

dfi \(k + ifi) l0 + k^ + i(M + iM) J 

A calculation along the line sketched above yields the simple result 

f d 3 k . . 

Pb = 2 I . 95 

Ja(k 2 )<fi (27r) 3 

This provides the familiar relation kj = 3n 2 pB verifying that this relation also holds 
in the context of the augmented model ([5]). 

We have studied the solutions M and M of the gap equations ( |91~1 , |92"|) as a 
function of the Fermi momentum kf. The result is shown in Figure [TT|. The coupling 
strengths Go and G c were chosen in order for the system to be in the confined 
phase (M ^ 0) at zero density. For increasing Fermi momentum kf one observes 
an increase in M up to a critical momentum kf, where M rapidly drops to zero 
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implying that the deconfined phase is realized at high density. The behavior of 
the color-singlet constituent quark mass M Q is different. It smoothly decreases and 
vanishes at kj. This behavior is different from that found at finite temperature (i.e. 
section |^) where the phase transition causes a discontinuity of Mq at the critical 
temperature. The deconfining phase transition at finite density is accompanied by 
the restoration of chiral symmetry as in the temperature case. 

Figure [12] shows the critical color-triplet coupling G^" 1 ^ as a function of the 
Fermi momentum kf. The dependence of G^ rit ^ on the Fermi momentum kf is 
qualitatively the same as the dependence on the temperature (compare Figure [ID). 
The critical density is nearly independent of the current mass of the quarks. 



7 Conclusions 



This paper describes how the NJL model can be modified minimally so as to 
take into account quark confinement. Confinement is understood in the sense that 
quark-anti-quark thresholds which plague the standard NJL model are screened by 
a random color background field: physical quantities are free of unwanted colored 
excitations. The mechanism that confines the quarks is analogous to Anderson 
localization in electronic systems and deconfinement can be induced by temperature 
and/or density in a way that seems to be consistent with QCD. In this model, 
deconfinement and chiral symmetry restoration occur at the same critical point. We 
have used, for simplicity, color SU(2) group but we have no reason to expect that 
qualitative features would be modified if we were to consider the realistic SU(3) 
gauge group. For applications to phenomenology, SU(3) gauge group will have to 
be treated. The extension to color SU{3) remains to be worked out. 

The problem that we hope to be able to resolve with the confining NJL model is 
to describe how hadrons - both mesons and baryons - behave in medium at finite 
temperature and density. Since confinement is suitably implemented in the model, 
we can treat the excitations that are not amenable to the conventional NJL model, 
such as scalar and vector mesons. For instance, we should be able to "derive" the 
BR scaling predicted in mean field of effective chiral Lagrangians P2| and address 
the properties of hadrons in hot and dense medium created in heavy-ion collisions 



or compact star matter p8| . A microscopic model of the kind presented here will 



have certain advantage over the "macroscopic" treatment of [32 
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Appendices 

A Some ingredients for the scalar correlation function 



Here we calculate the functions 11° (p 2 ), K a (p 2 ) and Il" /3 (p 2 ), defined in ([I6|jl7| , p!8| 



that enter in the scalar correlation function ([13]). For this purpose, we need to 
evaluate the trace, which extends over Lorentz- and color-space, of the functions 

tr {S{k + p)S{k)} , tv{r a S(k+p)T^S(k)} , tr {r a S {k + p) S {k)} , 

(A.l) 

where the fermion propagator S(k), (|T9"D, 

S{k) := f + i (MoV iM«t«) (A ' 2) 

possesses a non-trivial color structure. The Lorentz trace is straightforward. The 
color trace is most conveniently evaluated by introducing the eigenvectors |±) of 
the matrix M a r a which form a complete set in color space giving rise to a specific 
representation of the unit element, i.e. 

l = X»<i|. (A.3) 

i=± 



For instance, the color trace of the first expression in (|A.1|) is 



E tr L {(i\S(k + p)\l)(l\S(k)\i)} , (A.4) 

ii=± 



where tr/, indicates the trace over Lorentz indices only. Using 



S ®W = TTiuk±m & ' (A - 5) 
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one obtains for (|A.4j) 

tr L {s(k + p)s(k)} + (M -> -M) , (A.6) 

where 

S (A;) = — with .1 M • /.U . (A.7) 

It is now straightforward to calculate U®(p 2 ) in ([TED: 

The Lorentz trace can be easily performed. Introducing Feynman's parametrization, 
one has 



TO/' 2\ 



where Q = a(l — a)p 2 + A 2 . After the substitution q = k + ap, we obtain the final 
result 



= ^ - #o(p 2 ) , (A.ll) 

with H${p 2 ) defined in ([24]). Since the integral in flA.10 ) is divergent, we cut off the 
q- integral by a sharp 0(4) cutoff A. This regularization procedure is included into 
the definition of our low-energy effective quark theory We have ensured (see section 
^) that chiral symmetry is not violated by the cutoff procedure. 

For the evaluation of K a (p 2 ), we proceed along the same line as sketched above. 
We first perform the color trace in the last expression of ( |A.1| ): 

£ tr L {(i\r a S(k + p)\l)(l\S(k)\i)} = £ tr L {(z\T a \l)s(k + p)s(k)(l\z)} . 

i,l=± i,l=± 

(A.12) 

From fl22|) and the orthogonality of the eigenvectors |±), we find 

K a ( P 2 ) = -i— J ^tr L {s(k + p)s(k)} - (M--M) • (A.13) 

For the further calculation of K a , one might employ directly the above results to 
finally obtain 

2^ ,M a f 1 f d\ f q 2 -Q 



K a (p 2 ) = -U— / cia / — { — ^ — (M —> —M) , (A.14) 

yF J M Jo J (2vr) 4 \ [ g 2 + g2]2 v > v ; 
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and therefore K a (p 2 ) = ^-H v (p 2 ) with H v (p > defined in 

It remains to show that M a is an eigenvector of the polarization matrix 11"^. To 
do this, we first evaluate the color trace of 



tr jr°S(* + p)r«^pS(*-)} = £ tr L j (i|r°|(>»(i + p)s(k){l \r 



M 1 ' 



(A.15) 

Since |±) are eigenvectors to t^M^ , the expression (|A.15 ) is 

— tr L + p)s(fc)} + (M — — M) . (A.16) 

This expression was already calculated to obtain 11° (p 2 ). Since M@ is also an eigen- 
vector of the interaction matrix G a/3 (0), the final result is 

nf(p 2 )^ = (^r-Ho(p 2 )) (A.17) 



M \G r y ^ V M 
B The equations of motion 



In this section we derive an explicit expression for the equations of motion (|TTD 
and fll2|) . To perform the color trace of the quark propagator S(k), we introduce 
the eigenstates |±) of the color matrix M a r a (compare Appendix [AJ): 



tr5(Jfe) = tr L X^'I^W V) 



where trj, is the trace over Lorenz- indices only. The equation of motion (|11|) becomes 



(M - m) - i 



d 4 k ( #-iA 



+ (M — > — M) 



0. 



(B.2) 



(2tt) 4 [P + A 2 

where A was defined in ( |A.7|) . It is straightforward to evaluate (tB.2|) . One obtains 

d 4 k k 2 + M 2 + M 2 



i_(M -m) - 8M / 

(jTQ J 



(2tt) 4 A; 4 + 2(M 2 - M 2 )k 2 + (M 2 + M 2 f 



. (B.3) 



Setting M a = iM a , we derive the explicit form of the remaining equation of motion 
(|l2j) in a similar fashion. The color-triplet part of the quark propagator, tr S(k)r a , 
is 

M a 

tr L J2(i\ S(k)r a \i) = tr L $>(£;) (z| r a |*> = tv L s(k)— - (M -> -M) , (B.4) 



i=± 
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with s(k) from ( |A.7|) . If M a is an eigenstate of the interaction matrix, i.e. G^M 13 = 
G c M a , the equation of motion (|T^) reduces to 

A direct calculation yields 

—M - 8M y fc4 + 2 (M 2 - M 2 )fc 2 + (M 2 + M 2 ) 2 ' ( ^ 

Both integrals in (p.3j ) and (|B.6|) can be expressed with the help of the functions Iq 



and I v defined, respectively, in fl36[) and (|37|) , to obtain the desired result presented 



in (33) and 



C Analytic structure of the scalar correlation function 



As argued in section |3]2|, we believe that the approximations made to derive the 
scalar correlation function do not violate any fundamental requirements of quantum 
field theory. Nevertheless we observe a non-trivial analytic structure of the scalar 
correlation function which might make one suspect that some constraints may be 
incorrectly implemented. In this subsection, we show by an explicit calculation 
that one particular constraint, namely, the stability criterion, is indeed satisfied. Of 
course there are many more constraints. Here we shall offer one evidence that our 
model is compatible with a general axiom of quantum field theory. The investigation 
of further constraints which might provide insight into the possible analytic structure 
of Green's functions of a confining theory seems very interesting to us and will be 
relegated to a future work. 

For completeness, we rederive the stability criterion. Following the standard pro- 
cedure |Tl], the correlation function of a local operator j(x) is written in momentum 



space as 

G(q) = [d 4 xe^ (0\Tj(x)j(0)\0) , (C.l) 



where T is the time-ordering operator. Rewriting the time ordering, one has 

G(q) = J d 4 x <*» 9(x )(0\ [j(x),j(0)} |0> + J d A x e** (0\j(0)j(x)\0) (C.2) 
Inserting a complete set of eigenstates, the last expression becomes 

,A X el ( g+Pn )x ( |j(0)|n)(n|j(0)|0) (C.3) 
= ^(2vr) 4 5 4 (g + Pn )(0|j(0)|n)(n|j(0)|0). 
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In the lab frame (q° > 0), there is no contribution from the term ( |C.3|) to the 
correlator ( Clj) , if we are dealing with a stable theory, since there are no states with 



negative energy {E n = p„ > 0). We therefore observe that the correlation function 
vanishes for xq less than zero, i.e. 

(0|Tj(x)j(0)|0) = f(x)6(x ). (C.4) 

In the standard case of a constituent quark model, the stability criterion is satisfied 
for the following reason: the correlation function in momentum space is analytic in 
the upper q complex-plane and possesses cuts at the real axes due to quark anti- 
quark thresholds (see Figure [13] (a)). For x < one might calculate the Fourier 
transform by closing the path by a semi-circle in the upper half plane. Since there 
are no poles or cuts, we conclude that the Fourier transformation yields zero. 

In order to check the stability criterion in our model, we first study the 
analytic structure of the scalar correlation function fl39l) . To this aim, we have 
explicitly calculated the functions H / V (p 2 ) in (^|) and (|25|). The result is 

H (p 2 ) = h(p 2 ;M ,M) + h(p 2 ;M ,-M) , 

H v (p 2 ) = -i{h(p 2 ;M ,M)-h(p 2 ;M ,-M)} (C.5) 
4tt 2 A 2 1 (AA 2 \ 3/2 + 1 




(C.6) 



i(^ + i)(^ + i) ln V(^ + 1 ) + 1 



+ 1-1 



2 y^ 4A 2 +4A 2 | ^ y^ 4A 2 +4A~ 

where A = M + iM. We confine the momentum to |p| 2 < A 2 . In this case the only 
cuts which occurs in the upper half go complex-plane are shown in Figure |13|(b). 
Their orientation is given by the phase <fi of the complex number M + iM. In order 
to perform the Fourier transform to the coordinate space, we close the path in the 
upper half plane by the contour depicted in Figure [l3|(b). It seems that the Fourier 
transform does not yield zero as it would be required by the stability criterion, since 
there are now the contributions from the cut. Our main observation is, however, 
that these contributions exactly cancel. In order to see this, we first note that the 
correlation function A(g ) possesses reflection positivity, i.e. 

A(z*) = A*{z) , (C.7) 

as the functions H / v do. This can be shown either by a direct inspection of ( |C.5| ) 
or, which is more instructive, by tracing it back to the fact that we consider color 
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singlet correlation functions. In the latter case, the appropriate superposition of 
terms with M and — M provides ( JCJ.7| ) . Equation (U.7) immediately implies that 



the complex part of the integration over the contour in the half-plane vanishes, since 
this contour is chosen to be symmetric with respect to the imaginary axis. In order 
to prove that the contributions from the cuts cancel, it is sufficient to show that 
the contribution from either cut is purely imaginary. That this is indeed the case 
holds on very general grounds. Although this may be well-known to the readers, we 
nevertheless feel that we should give some arguments on this matter. We have to 
prove that integrals of the type 

I := J dxK(\n(x))f(x) (C.8) 

produce purely imaginary results, where C is the contour surrounding the cut (like 
the contour ACB in Figure |l3|(b)). In fact, one has to study the more general 
integrand K[ln(g(x))]f(x), the integral of which, however, can be traced back to 
the one in (|C.8|) by a change of variable x — > y = g(x) and a redefinition of the 
function f(x). This change of variables provides us with a linear cut in the complex 
plane. If we define the cut of the logarithm to coincide with the negative real-axis, 
a straightforward calculation yields 

/= f° dx fix) {K(\n(-x) + in) - K(\n(-x) - in)} , (C.9) 

J-R 

where R is the length of the cut under considerations. To continue we have to 
assume that the function K(x) possesses a Laurant-expansion around x — 0. It is 
straightforward to check that this is the case in the context of the scalar correlation 
function of our model. We then have 

oo o 

/= Y, c n dxf{x) {(ln(-x) + m) n - (ln(-x) - in) n } . (CIO) 

n=-m , '^ R 

The even powers of in drop out and we end up with a purely imaginary result. This 
completes the proof. 

To summarize our results, the scalar correlation function possesses a nontrivial 
analytic structure in momentum space, e.g. consisting of cuts in the upper half q 
complex-plane. We have proven that the particular analytic structure is compatible 
with the stability criterion (|C.4|) due to an intrinsic cancellation mechanism of the 



contributions from the cuts. 

We finally present the scalar correlation function in Minkowski-space. From 
the very beginning, our model is given in Euclidean space and the Green's functions 
in Minkowski-space are defined by the Wick rotation (see Figure [T3|). In our case, 
this Wick-rotation is non-trivial, since we have to take into account the contribution 
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from the cuts in Figure [14]. Suppose the scalar correlation function in Euclidean 
space is part of some scattering amplitude, e.g. 



/oo 
dtfi A(q 2 )g(q 2 ,p 2 ) (C.ll) 
-co 



where Q{q 2 ,p 2 ) is assumed to be analytic in the first and third quadrants of the 
g,f -plane. From Figure [14| we have 



/ioo roo 
dq* A(q 2 ) g(q 2 ,p 2 ) + 2i / dx I cut , (C.12) 
-ioo J—oo 

where we have used that the contribution I cut from one cut is of the form ( |C.10| ) 
and purely imaginary (I cut is real). Performing the substitution q^ 1 = —iq^ and 
redefining x = q(f , we obtain 

/oo roo 
dq™ A{-q 2 ) G(-q 2 ,p 2 ) + 2i / dq^ 1 I cut , (C.13) 
-co J—oo 

and arrive at the scalar correlation function in Minkowski space 

A M (q 2 ) = A(-q 2 ) + 2I cut (C.14) 

Note that the contribution from the cuts is real. We expect I cut to contribute only 
a background to A M (q 2 ), whereas the significant structure, such as imaginary parts 
from thresholds and poles from particle states, is produced by the Euclidean scalar 
correlation function at negative momentum squared. 

D Reduction of the Bethe-Salpeter amplitude 

In order to solve the Bethe-Salpeter equation ([|1]) for the hidden color structure 
of the pion, one has to perform the color trace of polarizations containing two 
quark propagators S(k). Since the calculational technique parallels very much that 
performed to obtain the scalar correlation function, we only sketch the derivation 
briefly and refer the reader to Appendix [A] for further details. The traces of interest 
are 

tr{ l5 S(k + p) lb S(k)} = tr L { l5 s(k + p) l5 s(k)} + (M - -M) , 
ti{r a l5 S(k + p) l5 S(k)} = ^-tr L { l5 s(k +p) l5 s(k)} - (M -> -M) , 

tr{T a l5 S(k + p)T? l5 S(k)}^- = ^ t T L { l5 s(k + p) l5 s(k)} + (M^-M) 

(D.l) 
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with s(k) defined in ( A..7). The next step is to calculate the trace over Lorentz 
indices tr^: 

k 2 + A 2 + kp 



# + i>-iA U-iA 



(D.2) 



' {k + p) 2 + A 2 10 k 2 + A 2 j [{k + p) 2 + A 2 } {k 2 + A 2 ) 

To obtain the desired matrix elements entering into (41), an integration over the 
loop momentum k is required. Introducing Feynman's parametrization and shifting 
the momentum integration q = k + ap yields 

Tr { l5 s(k + ph 5 s(k)} = -4 f da [ (q) - { ~ ° { } - + 2 . (D.3) 



[q 2 + a(l - a)p 2 + A 2 } 2 
Inserting ( p.3|) into ( p.l[ - P.lD , we have all the ingredients to derive the result (fE| 



E Calculation of the electromagnetic form factor 



Here we evaluate the matrix element (|57"D, which is directly related to the elec- 
tromagnetic form factor. For this, we first perform the color trace by the method 
employed several times before (for details see Appendix Q): 



{p 2 - pD wl - 2P.P, w; . 



(E.l) 



The functions 



-M) 



(E.2) 



W v 



' I 7^7 tCL (75 S(k - I) 7 , S(k + |) 75 S(k - P - ) \ - ( \ I \ f ) 



(2tt)< 



(E.3) 



are introduced for abbreviation. The propagator s(k) is given by ( |A.7|) . Performing 
the Lorentz trace of the terms under investigation in ( |F.2| , |F.3D , i.e. 

+ P -f ) 



(k - |) 2 + A 2 (k + |) 2 + A 2 



(E.4) 



one obtains 



w v 



-Pu / ^7 t i-f T + (M —* -M) (E.5) 

2^ J (2tt) 4 (A; _£)2 + A 2 (fc + f) 2 + A 2 1 > \ ) 



2> J |_ V 2/ 

i r d 4 k 1 



(2tt) 4 (/c - |) 2 + A 2 (fc + f) 2 + A 2 



(M -> — M) (E.6) 



34 



It is of a particular interest to expand these functions for small p 2 because they are 
required to normalize the pion Bethe-Salpeter amplitude and to calculate the pion 
charge radius. A direct calculation yields 



W 0/v 



(F 0/v - P 2 Ro /v + 0{p 4 



where F Q / V are defined in ( f47[ , f48D and 



R 

R v 



d A k A 2 - k 2 



(2tt) 4 (k 2 + A 2 ) 
d A k A 2 - k 2 



(2tt) 4 (k 2 + A 2 ) 



) + (M — > —M) , 
- (M -> —M) 



(E.7) 

(E.8) 
(E.9) 



Inserting ( |E.7p in ( |E.1[ ) leads to the desired matrix element. 

In order to obtain the electromagnetic form factor for non-vanishing momentum 
in Minkowski space, an analytic continuation of the Euclidean momentum squared 
to negative values is needed. This analytic continuation is difficult from a technical 
point of view, since the momentum integration in (|E.5| , |E.6| ) can be performed only 
numerically. Introducing 



K ± = 

a term of interest is 



(*±|) 2 



p2 

k 2 — ± kp 



A, 



M ±iM 



1 



[K- + A\] [K+ + A\] 



+ {A, 



A. 



(E.10) 



(E.11) 



We first remove the complex parts that enter via A±, namely, 

K + K^ + B_(K- + K+) + B 2 _- AM 2 M 2 
[K 2 _ + 2K_B_ + Bj] [K 2 + 2K + B_ + B 



21 ' 



(E.12) 



where B± = Mq ±M 2 . Introducing the angle a between the Euclidean four- vectors 
k and p, the crucial observation is that in ( E.12|) , only the terms quadratic in the 
external momentum p appear: 



K+K_ 



Kl + K 2 



K + + K. 



Il l P \2 i 2 2 2 

(k — —J — k p cos a 



(k 2 



T 

2(k 2 - P - 



,7 2 2 'A 

+ k p cos a 



(E.13) 
(E.14) 

(E.15) 
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It is now a simple matter to perform the analytic continuation p 2 — > — p 2 . We finally 
present the result for the functions W°, in ( |E.5| , |E.6| ) 



1_ /• g_ + 2T_/3 + T 2 - AM 2 M 2 

^3 P/t 7 a 2 + AT_a.f3 + 2T 2 a + + AT 2 a. + Aj 

L p 4M M(/3 + T. 

" «7r 3 ^ J a 2 _+ AT_aB + 2T 2 a+ + AT 2 c- ^ ^ ^ 2 * ^ ^ 4 



M ^ :5 ''" / '" ™ 2 -l 4T_a_/3 + 2T?a + + 4T 2 ct_ + AT_Tlf3 + Tf ' 



8vr 3 ^ 7 a 2 . + 4T_a_/3 + 2T 2 a + + AT 2 a_ + AT_T 2 f3 + T\ 

(E.16) 



where dV = da sin 2 a d{k 2 ) k 2 , a G [0, 7r] and 

2 2 

a± = (k 2 ± ^-) 2 ± k 2 p 2 cos 2 a , /3 = k 2 - V — , T± = M 2 ± M 2 . (E.17) 

The momentum integration k 2 and the angle integral in ( |E.16| , [ETT6|) are left to a 
numerical calculation. 
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Figure Captions 



Figure 1: The color-singlet and color-triplet constituent quark masses M and M 
for H^t- = 2 and m/A = 0.01 for the confining phase (solid and short dashed line) 
and for the deconfined phase with M = (long dashed line). 



Figure 2: The critical color-triplet coupling G c as function of the color-singlet 
coupling G for different values of the current mass m in units of the cutoff A. 



Figure 3: The scalar correlation function as function of the Euclidean momentum 
transfer for different values of the current quark mass m for ^4- = 1 and ^f^- = 2. 
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Figure 4: Graphical representation of the electromagnetic form factor (a) and the 
pion decay constant (b). 

Figure 5: The electromagnetic form factor F{p 2 ) as function of the negative Eu- 
clidean momentum — p 2 . 

Figure 6: Mass m n and decay constant f n of the pion as function of the current 
mass m for = 1 and = 2. Also shown are the constituent quark masses 
M and M. 

Figure 7: The pion dispersion relation as function of the Euclidean momentum 
transfer for zero current quark mass m = and for = 1 and = 2. 

Figure 8: The pseudoscalar correlation function as function of the Euclidean mo- 
mentum transfer for a current quark mass m = 0.01A and for = 1 and ^£0- = 2. 

Figure 9: The color-singlet and color-triplet constituent mass as function of the 
temperature T. M , M, and T in units of the cutoff A. 

Figure 10: The color-triplet critical coupling strength G c as function of the tem- 
perature T for different values of the current mass m in units of the cutoff. 

Figure 11: The constituent quark masses Mq and M as a function of the Fermi 
momentum kf for = 1 and = 2. M , M and kf in units of the cutoff A. 

Figure 12: The color-triplet critical coupling strength G c as a function of the Fermi 
momentum kf for different values of the current mass m in units of the cutoff. 

Figure 13: The analytic structure of the Euclidean scalar correlation function in 
the complex q plane for the standard NJL model (a) and the confining model (b). 
Cuts are indicated by dashed lines. 

Figure 14: The Wick rotation to relate the scalar correlation function in Euclidean 
space to that in Minkowski space. 
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